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Abstract 



This article considers the semilinear boundary value problem given by 
the Poisson equation, — Am = f(u) in a bounded domain C R n with 
smooth boundary. For the zero boundary value case, we approximate a 
solution using the Newton-imbedding procedure. With the assumptions 
that /, /', and /" are bounded functions on R, with /' < 0, and fl C R 3 , 
the Newton-imbedding procedure yields a continuous solution. This study 
is in response to an independent work which applies the same procedure, 
but assuming that /' maps the Sobolev space H 1 (fi) to the space of Holder 
continuous functions C°"(Q), and f(u), f'(u), and f"(u) have uniform 
bounds. In the first part of this article, we prove that these assumptions 
force / to be a constant function. In the remainder of the article, we prove 
the existence, uniqueness, and if 2 -regularity in the linear elliptic problem 
given by each iteration of Newton's method. We then use the regularity 
estimate to achieve convergence. 

Introduction 

The goal of this article is to find suitable hypotheses on a function / S C 2 (M.) 
related to attaining a solution to the semilinear boundary value problem given 



using the Newton-imbedding procedure that is applied in [2]. Here, /(it) is 
defined as / o u. In this sense f can be viewed as a map from a space of real- 
valued functions to another space of real valued functions via composition. In 
addition, H k (il) is defined as the L 2 functions on having (weak) i th derivatives 
(1 < \i\ < k) which are L 2 functions on f2. This is the Hilbert space notation 
substituted for the Sobolev space notation W k > 2 (n). The space of real- valued 
functions on f2 which are Holder continuous with exponent a will be denoted 
C Q (f2). The author of [2 J achieves an H 2 solution when fHs a domain in M 3 
and T is smooth, provided the following assumptions on / hold: 

• 1. / is a continuous map from H 2 (£l) to L 2 (Q). 

• 2. /' and /" are continuous maps from iJ 1 (SI) to C Q (f2), a <G (0, |]. 



by 




-An - f(u) 



in 17 

on r = on, 



• 3. There exists a constant M > such that 



\\f(u)\\mn) < M for all u £ H 2 (n), \\f(u)\\ cam < M for all u £ H^Sl), 



• 4. (— /') is positive in the sense that (—f'(u)v, v) > for all =/= v £ H 2 {£1) 

An additional condition in [5] is the choice of a uniform width of time intervals 
in the procedure that ensures convergence, which exists as a consequence of the 
above assumptions. However, we prove the following theorems in Sections 2 and 
3 of this article: 

Theorem. 2.1 If f : R — » M is a map from to C°(Cl) via composition 

and f2 is a domain in R n with n > 2, then f is a constant function. 

Theorem. 3.1 Let h : K -> K map i7 2 (fi) f| #o(fi) to Z,P(fi) via composition, 
where 1 < p < oo and is domain in R n . If there exists a constant M > such 
that \ \h(u)\\LP < M for all u £ H 2 (fl) f]H^(fl), then h is a bounded function 
on R, i.e. there exists a constant C > such that \h(x)\ < C for all ieR. 

By Theorem 2.1, the assumption in (2) that /' maps H 1 to C Q forces /' to 
be a constant function. Theorem 3.1 shows that the uniform bound on f(u) in 
assumption (3) forces / to be a bounded function on R. Thus / is shown to 
be linear and bounded on R, and is therefore a constant function, reducing the 
scope of the procedure in [2 to the family of problems given by — Au = const. 

In Section 1 of this article, we construct a 'mesa' function (see Figure 1 in 
Section 1) whose existence in i/ 1 (SI) will serve as a counterexample to a non- 
constant mapping. In Section 2, the mesa function is used to prove Theorem 
2.1. In Section 3, Theorem 3.1 is proven using a sequence of smooth 'bump' 
functions in H 2 . As a consequence of this, the uniform bounds also imposed 
in (3) on f'(u) and f"(u) imply that /' and /" are also bounded functions on 
R. In Section 4, we describe and apply the Newton-imbedding procedure to 
the case of (*) with a zero boundary condition. Of primary importance in the 
procedure is the following linear boundary value problem, 



given by each iteration in the Newton-imbedding procedure. Here, q(x) is a 
positive scaling of (— /') while g(x) depends on / and /' in a manner that 
allows g £ L 2 under our assumptions. The exact hypotheses on q and g will 
be made precise in Section 4. As in the assumption that q > allows for 
existence and uniqueness for (**) in H 1 , as well as the regularity lifting of the 
H 1 solution to H 2 . For the remainder of the article, it will be understood that 
(**) is the general boundary value problem stated above, with the conditions 
that g £ L 2 and q > 0. Under the following assumptions on /, 



and ||/"(u)||c°(o) < M for all u £ H 1 ^). 




on r, 



in S7 
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• I. / is a continuous map from H 2 (Cl) to L 2 (Q). 

• II. /' and /" are continuous maps from to L n (fl) 

• III. there exists a constant M > such that 

|/|<M, |/'|<M, and |/"| < M. 

. IV. (-/') > 0, 

we prove existence and uniqueness for (**) in Section 5, and achieve the reg- 
ularity lifting of an Hq solution of (**) to H 2 in Section 6. These results are 
summarized in the following theorem: 

Theorem. 6.1 Let VI be a bounded domain in M. n with smooth boundary T 
and n > 2. Then for g e L 2 (fl), q e L n (Q), and q > 0, the linear boundary 
value problem 

, , J —Au + q(x)u = g(x) in ft 
^*'{ u\ r = onT, 

has a unique solution u e H 2 (fl) f] Hq with 

\\u\\m(Q) < C{\\g\\ L 2 (n) ), 
where C depends only on fl, n, and q. 

In Section 7, under an additional assumption (V) concerning the uniform 
width of time intervals in the procedure, convergence in the procedure is achieved 
resulting in the following theorem: 

Theorem. 7.1 With ft a bounded domain in M 3 with smooth boundary and 
assumptions (I)-(V), the semilinear boundary value problem, 

, f — Au = f(u) in 

u\ r = onT = dn, 

has a unique solution in H 2 (fl) f] Hq (Ct), and hence a continuous solution, 
which can be approximated by the Newton-imbedding procedure. 



1 The Mesa Function 

Let be a domain in R n with n > 2 and let c € CI. Since the function will be 
radially symmetric about c, define r =\ x — c | for a; e f2, and T > such that 
B(c,T) CC ft, where B(c,T) denotes the open ball of radius T about c. Also 
let a,b e R with a < b, and a G (0, n ^ L ). In order to define the function, it is 
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necessary to decompose the interval [0,T] as follows: 



If we let = —, then there is an such that 



1 



1 



(«+)" (r+)« 



b — a. 



In particular, < < r±. Setting s 1 = allows for an r x such that 



1 



(rr) Q (sr) c 



= b — a. 



In particular, < r 1 < s x . Continuing in this manner, set r+ +1 
Note that rlt , -, > for all m and r\t , , goes to zero with 



Using the above notation, let : f2 — > R be the radially symmetric piecewise 
function defined inductively by 



r 



U(r) = < 



+ a 



r>T 

r+ < r < T 

s+ < r < r+ 

s m ^ r — s m 

r~ <r<s~ 



' m+1 



<r < r . 



We will call U(r) a mesa function with exponent a. Figure 1, below, is a 
sketch of a mesa function whose partition points have been altered to show more 
'mesas'. 



U is bounded and has compact support, so is trivially in L 2 (fl). It remains 
to show that it has (weak) first derivatives in L 2 (fl). The proposed first deriva- 
tives are given by 
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r 

Figure 1: Artist's depiction of a mesa function 



U Xi (r)={ 






, r > T 


-2a 
T 


, r+ <r <T 


— axi 

r a + 2 


! s m ^ r ^ s m 





i s m < r — s m 


axi 

r a+2 


, r" <r<s m 





, r+ +1 <r<r 



Away from zero, on each annulus of the decomposed fl, the expressions in 
U Xi are classical derivatives of their corresponding expressions in U(r). Let 
4> G Co°(f2) and fix AT. Integrating U<p Xi by parts over the annuli given by 
[rt-T], [s£,r+], [ s m, s mL [ r m> s m]> and [^m+i^m] for m = h -, N and recall- 
ing that {/ = for r > T, gives 

/ U(f> Xz dx=- I U Xi (j)dx+ / U^dS, 

Jn-B(c.r+ + 1 ) Jn-B(c,r+ +1 ) JdB(c,r+ + 1 ) 

where p = (p 1 , ...,p n ) is the inward pointing normal on dB(c, r^ +1 ). 
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Let u(r) = —. Note that \U Xi \ < \u x .\, so that \DU\ < \Du\. 

Following the line of argument [1, p. 246] given by L. Evans, since a < n — 1, 
\Du\ = G L^fi) and therefore |£>[/| G L^fi). 

Letting N ^ oo (and thus r^ +1 — » 0), 

< || Ety Hoc / < M(r+ +1 )"" 1 - 0, 

Jas( c ,r+ +1 ) 



hence 




Therefore t/^ is a (weak) derivative of U. Moreover, since a < — - — , 

following the argucmcnt in [1, p.246], | Du |g L 2 (tt) and thus \DU\ e L 2 (tt) 
and U(r) G _ff 1 (fi). The following lemma summarizes the above discussion: 

Lemma 1.1 If ft is a domain in M. n with n > 2, and U(r) is a mesa function 
n — 2 

with exponent a < — - — , then U(r) G _ff 1 (fi). 



9B(c,r+ + 1 ) 



2 Constant Mapping 

Theorem. 2.1 If / : R — > M is a map from ii^ 1 (SI) to C°(f2) via composition 
and £1 is a domain in R™ with n > 2, then / is a constant function. 

Proof. Suppose on the contrary, that / is not constant and assumes distinct 
values at a and b. Without loss of generality, assume that a < b. Let ceSJ and T 
be such that B(c, T) CC ft. Since n > 2, there exists a such that < a < r± ^ L . 
Let U(r) be the mesa function centered at c, with exponent a, support in 
B(c,T), and prescribed maximum and minimum, b and a, respectively. By the 
above lemma, U(r) is in _ff 1 (fi). Using the notation in the previous section for 
the domain of U(r), it holds that for any 5 > there exists an N such that 
[s w ,s+] C B(c,6) and [r^ +1 ,r JV ] C B(c,S). Note that foU = f(b) on [s w ,s+] 
and / o U = f(a) on [r~^ +1 ,rjf]. Since the measure of the above intervals is 
strictly positive, / of/ has no continuous representative. In other words, the 
oscillations of / o U do not diminish in any neighborhood of c. This contradicts 
the hypothesis that / maps U to a continuous function. 

□ 

Now, as an immediate application of Theorem 2.1, the assumption in (2) 
that / maps H 1 into continuous functions forces / to be constant. 



G 



3 Uniform Bounds 



For this Section we assume £1 is a domain in K". 

Theorem. 3.1 Let f : R R map H 2 (n)f]H^(n) to LP (SI) where 1 < 
p < oo. If there exists a constant M > such that \\f(u)\\i,p < M for ah 1 
it G H 2 (fl) P| Z/o(fi), then / is a bounded function on R, i.e. there exists a 
constant C > such that |/(x)| < C for ah x £ R. 

Proof. Let p < oo. Suppose on the contrary, that / is not bounded. Then there 
exists a sequence, {a^}^! in R such that |/(£fc)| > k. Let y n E il and r such 
that B = B(y ,r) CC ft. Set i?i = -B(yo> §)• Choose a smooth function, 7, 
such that 7 = 1 on Bi, 7 = on il — B, and < 7 < 1. Define the smooth 
function Uf. on by = Xkj- Then Uk € H 2 (il) f] Hq(£1) for all fc and 

||/K)|U P( n) > 11/(^)11^(50 - 11/(^)11^(50 > k\Bi\5. 

2 2 z 

Choosing fc large enough such that fco|-Bi| ? > M gives a contradiction. If 
p = 00, a similar computation holds, choosing fco > M. 

□ 

Remark: Since the C a norm has the L°° norm as a summand, Theorem 5.1 
with p = 00 suffices to show that a uniform bound on ||/(u)||c° implies / is 
bounded. Therefore the assumptions made in [2], imply that /, /', and /" are 
bounded functions. Moreover, under the same assumptions, as shown in the 
previous Section, / is linear. In this case / is a constant, reducing the scope of 
the procedure to problems given by — Ait = const. 

4 Newt on- imbedding Procedure 

The Newton-imbedding procedure we wish to apply to 

{ is / -Aw = f(u) in tt 
1 } \ u\ r = onT = dn, 

has two parts. It is well described in [2j, but recalled here for clarity. The 
procedure first imbeds the problem in a one-parameter family of problems, 

—Am = tf(u) in il 
with u = on L and parameter t £ [0,1]. We set 

F t (u) = Au + */(«). 

Solving (*') is then a matter of solving Fi(u) — 0. Let u(x, t) be the solution to 
F t (u) — 0. Starting with to = 0, the problem is solved with solution u(x, 0) in f2. 
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Observe that with boundary value zero imposed, u(x,0) is uniquely determined 
as u(x,0) = 0. To solve F tl (u) — 0, u(x,0) is taken as an initial approximation 
and the standard Newton's method is applied. With convergence, the solution 
u(x,ti) to F tl (u) — is achieved. The function u(x,ti) is then used as an 
initial approximation for F t2 (u) = and so on for increasing times tj . Thus the 
solutions are pushed along with increasing times using Newton's method with 
the goal of reaching t — 1 in finitely many time shifts. Let uo(x, tj) — u(x, tj—i), 
the initial approximation for F tj (u) — and u m (x, tj) be the m th iteration of 
Newton's method at time tj. In the following discussion, the argument of the 
u m 's will be suppressed. We will also temporarily use the symbol D for the 
Frechet derivative in contrast to its usual use as the gradient. Note that 

DF tj (u m )[w] =Aw + tjDf(u m )[w] and Df(u m )[w] = f(u m )w 

for w € H 2 (Q) and that the (m + l)*' 1 iterate in the Newton approximation is 
given by 

DF t .(u m )[u m+1 - u m ] = -F t .(u m ). 

In this case, the (to+1)*' 1 iteration at time tj yields the following linear problem: 

/ f -Au m+ i + (-tjf'(u m ))(u m+ i) = tj(f(u m ) - f'(u m )u m ) in 
1 ! \ u m+1 \ r = onT. 

This is the problem 

-Au + q(x)u = g(x) in 
{ >\ «|r = onT, 

stated in the introduction with 

1 = -tjf ( u m), 9 = tj[f(u m ) + f (u m )u m ], and v = u m+l . 

Initially, a weak solution in Hq is desired, so it makes sense that u be in 
Hq and that / and /' should be defined on Hq. However, as will be shown in 
Section 6, an Hq solution to (**) is also in H 2 . In light of this, / and /' need 
only be defined on H 2 . Note that if / maps H 2 to L 2 and /' maps H 2 to L n , 
then g is in L 2 for all dimensions n > 2, via the Sobolcv imbedding theorem. 

1 271 

Indeed, since u is in iJ , u is again in L™- 2 and the Holder inequality gives 

/ lf'(u)u} 2 <C\\f(u)\\U\u\\ 2 ^ 
Jn Ln ~ 2 

To fullfill the positivity condition on q in (**), we impose that — /' > 0. Now, 
at each time tj > and for all m, the m th step in the iteration at time tj is a 
model for (**). 

For the remainder of the article, we assume is a bounded domain in M. n>2 
with smooth boundary T and make the following assumptions (I)- (IV) on the 
nonlinear function /: 
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• I. / is a continuous map from H 2 (£l) to L 2 (il). 

• II. /' and /" are continuous maps from to 

• III. there exists a constant M > such that 



\f\<M, \f'\<M, and |/"| < M. 



. IV. (-/') > 0. 



Remark: There is a redundancy and lack of 'sharpness' in assumptions (I) 
and (II), given (III). Indeed, if the functions /, /', and /" are bounded, they 
naturally map to bounded functions on f2, and hence to L°°(f2) which is con- 
tained in L p (il) for all p > 1 since £1 is bounded. The reason for stating L 2 
explicitly is that it is a familiar assumption for framing weak solutions to linear 
elliptic problems. The bounds on the functions are not necessary to existence 
and uniqueness in (**), nor to the regularity lifting of the solution to H 2 
. Moreover, the L 2 hypothesis on / and the L n hypothesis on /' are suffi- 
cient for existence and uniqueness and the regularity lifting. For a more general 
treatment of elliptic equations with measurable coefficients, see [3]. 

5 Existence and Uniqueness 

For this Section, we assume (I), (II), and (IV). To prove existance and unique- 
ness for (**) in Hq(£1) (H 1 functions with zero on the boundary), the Riesz 
Representation theorem is sufficient. We seek a unique solution in Hq(Q). The 
associated energy form for (**) is 



It is well defined on Hq(H). Indeed, since n > 2 and u,v S Hq(Q), then 
it, v G L"^ 3 (f2) by the Sobolev imbedding theorem. Also since O is bounded, if 
q £ L n (ft), then q e Lt(fi). Note that 




£2 



DuDv + quv. 



2 



n-2 



n-2 



+ 



+ 



= 1. 



n 



2n 



2n 



Therefore by Holder's inequality, quv is integrable over Q with 




This inequality combined with the Sobolev inequality 




<C\\u\\ H i 



gives 
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\B(u,v)\<C\\u\\ H i\\v\\ Hi 

where C > is dependent on f2, n, and ||<z||r§ but not on u and v. By the 
Poincare inequality and the positivity of q, we have 

\\u\\ 2 H i<C I \Du\ 2 <C I \Du\ 2 + qu 2 = CB(u,u) 
Jn Jn 

where C > is dependent on n and f2 but not on u. Since / S L 2 (f2), it is a 
bounded linear functional on 77g(f2) [P. Since B(u,v) is an inner product on 
Hq, the Riesz Representation theorem provides a unique it* € Hq(H) such that 

= / fv for all u € #o- 

In other words, u* is the unique weak solution to (**) in Hq. 

6 Regularity 

With the same hypotheses as in the previous Section, we wish to lift the regular- 
ity of the unique solution to (**) from Hq to H 2 , with the estimate controlled 
by the L 2 norm of g(x). Theorem 6.3.4 (Boundary £f 2 -regularity) in [T] gives 
the desired regularity lifting of a solution to (**) when q 6 L°° . However, the 
L°° condition is only used in factoring out from the following integral to 

find, for u, v £ H 1 and e > in Cauchy's inequality, 




The L n hypothesis on q provides, 




by Holder's inequality, the Sobolev imbedding theorem and Poincare's inequal- 
ity. By the above estimates, we have also 

J {cuf < M\\Du\\ 2 L2 . 

Following the line of reasoning in [1. , the result for q € L n is a sufficient replace- 
ment for the estimate for L°° to get the regularity estimate, 

IM|jj2 (n ) < C (||s||L2( n ) + ||u||fl-i(n)) 
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where C depends only on f2 and n and q. Now, recalling the second energy 
estimate above, 



lul&Wfi) ^ = C Jju<C Q|b||| 2(Q) + ^|M| 2 



L 2 (fl) 



since u is a weak solution to (**). The last inequality is given by Cauchy's 
inequality with e = 1. Also since u is a unique solution, the L 2 norm of u is 
controlled by the L 2 norm of g by Theorem 6.2.6 in pQ. Therefore, 

IMI// 2 (0) < C\\g\\ L 2 (n) , 
where C depends only on fi, n, and more significantly, q. 

To summarize the results in Sections 5 and 6, we have: 

Theorem. 6.1 Let be a bounded domain in E™ with smooth boundary T 
and n > 2. Then for g £ L 2 (il), q e L n (f2), and g > 0, the linear boundary 
value problem 

, > J — Au + q(x)u — g(x) in f2 
\ u| r = on T, 

has a unique solution u £ H 2 (il) f] Hq with 

\\u\\Hi m <C(\\g\\ LHn) ), 
where C depends only on f2, n, and g. 

7 Convergence 

In the previous two Sections, it was shown that (**) is uniquely solvable in H 1 
and the solution is a priori in H 2 with estimate controlled by the forcing term 
g. Recalling that (**) represents an arbitrary iteration of Newton's method at 
time tj, the linear equation solved by the difference, u m +i — u m for m > 1, is 
given by 

-A(u m+ i - Um) + (-tjf '(u m ))(u m +i - u m ) 
= tj(f(u m ) - f(u m -i) - f(u m -i)(u m - u m -i)) in n 
u m +i -u m = on T. 

This is (**) with 

V = U m +l - U m , q=-tjf'(u m ), 
g = tj(f(u m ) - f(u m -l) - f'(u m -l)(u m ~ "m-l)), 
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and a zero boundary condition. Indeed, using the same argument as at the end 
of Section 3, it is clear that j e L 2 . For m = 0, by the definition of uq at time 
tj, the problem satisfied by u\ — uq is 

-A(ttl - M ) + {-t 3 f(u Q ))(ui - U ) 

= (tj -tj-i)f(u ) in fi 
ui — uq = on r, 

and is again a model for (**). To facilitate the convergence estimates to follow, 
it will be helpful to use Taylor's theorem to simplify g. Similar to the application 
of a mean value theorem used in [S], for m > 1, g can be written as 

9 = tj(u m - u m _i) 2 / f"(Tu m + (1 - T)u m _i)(l - t)cZt. 

J {0,1) 

Theorem 6.1 and the boundedness of /" give the estimate, 

\\u m +i-u m \\ H 2 < C\\tj{u m -u m ^i) 2 / f"{TU m + (1 - r)u m _i)(l - T)dr\\ L 2 

J (o,i) 



ct j . . , . o i 

< — i — \\{u m - u m _i) | 



L 2 



Ctj iVf ... \ 1 1 2 

< — Y~ \\{u>m-U m -l)\\ L *- 

Before progressing with the estimate, it is important to discuss the depen- 
dence on dimension. For dimensions n = 3 and n = 4, the L 4 norm is controlled 
by the H 1 norm, by the Sobolcv imbedding theorem, which in turn is controlled 
by the H 2 norm. For dimensions n — 5, 6, 7, and 8, the L norm is controlled by 
the H 2 norm, via the more general Sobolev inequality [l,p.270]. The subsequent 
calculations do not depend on which dimension n 6 (3,4,5,6,7,8) is assumed. 
However, only in dimension n — 3 does the general Sobolev theorem assure that 
our H 2 solution is indeed continuous. For n = 5, 6, 7, and 8, the H 2 solution is 
respectively, L la , L 6 , Lt, and L A . To continue with the convergence estimate, 
for n G (3,4,5,6,7,8), we have 

Ct j ]\^[ ... . 1 1 2 Ct j C s ... , . . 2 

1 IK-Wm - U m _i)|| i4 < \\{U m - U m _i)|| H a 

where C s is the constant from the Sobolev theorem and only depends on and 
n. Since in Theorem 6.1, C depends on \ \f {u m {x,tj))\\i,n and hence m and tj, 
we invoke the boundedness of /'. Therefore ||/'(it m (a;,fj))||.L» is bounded by 
some constant C > 0, uniformly over m and tj. Let K — CA ^ C " , Inductively, 

||u m+ i - u rn \\ H 2 < (tjK\\ui - u \\ H 2) 2 ~ 1 \\ui - u \\ H 2 



12 



and therefore for s G N, 

\\u m+s - u m \\ H 2 < [a 2 + ~ 1 + ... + a 2 - Uq\\ H 2 

where a = ijifHui — uollff 2 - If tj is chosen such that a < 1, then the positive 
expression in brackets above is bounded from above by the tail end of a con- 
vergent geometric series, and therefore goes to zero as m — > oo. We have now 
shown that Urn IS cl Cauchy sequence in the Banach space H 2 (£l), and therefore 
converges to some u* £ H 2 (il). As stated in [2], due to the continuity of / and 
the boundedness of /', it is clear that u* satisfies 

, ,s J — Au = tjf(u) in il 
^ ' \ u\ T = on T 

almost everywhere and that the uniqueness of the solution u* follows from the 
uniquenss of the solution u rn (x,tj) to (**) for each m and tj. One additional 
assumption is necessary for tj to be chosen as above, as well as for progressing 
to t = 1 in finitely many applications of Newton's method. Assumption (V) will 
be a condition on the width of the time intervals tj — tj-\. To make this precise 
we look at the the problem satisfied by U\ — Uq at time tj and apply Theorem 
6.1 and the boundedness of / and /' to estimate, 

IK - U \\ H 2 < C\\(tj - tj-l)f(u )\\ L 2 

< C(tj - i^OH/MH^ < MC(tj - tj-x). 

If A = MC, then A depends on the bounds on / and /', the volume of f2, and 
n, but not on tj. In the following inequality, 

Ktj\\u\ — uq\\h2 < KAtj{tj ~ tj-i) < 1, 

the condition for convergence was that the leftmost expression be < 1. Since 
tj < 1 for all j, it suffices to make the assumption (V): 

• V. For each j > 1, tj — tj-i < -j^j 

As K A only depends on il, p — 2, n, and M, (and in particular, not j), KA 
gives a uniform bound on the time intervals, and therefore t — 1 is attainable 
after finitely many applications of Newton's method. When is a domain in M 3 , 
the H 2 solution is then continuous by the general Sobolev imbedding theorem. 
We now list assumptions (I)-(V) and state the main result. 

• I. / is a continuous map from H 2 (fl) to L 2 (fl). 

• II. /' and /" are continuous maps from H 1 ^) to L n (il) 

• III. there exists a constant M > such that 

\f\<M, \f'\<M, and |/"| < M. 
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. IV. (-/') > 0, 

• V. For each j > 1, tj — tj-i < 

Theorem. 7.1 With ft a bounded domain in R 3 with smooth boundary and 
assumptions (I)-(V), the semilinear boundary value problem, 

, ,J — Au = /(it) in il 

u\ r = onT^dn, 

has a unique solution in H 2 (il) f] Hq (fi), and hence a continuous solution, 
which can be approximated by the Newton-imbedding method. 



8 Conclusion 

The goal for improving this procedure is to weaken the assumptions on / and 
/'. In particular, to eliminate the boundedness or equivalently the uniform 
boundedness of f(u) and f'(u). To do this requires a function / such that 
f(u m (x,t)) does not grow too fast in L 2 norm as t increases and such that 
f'(u m (x,t)) does not grow too fast in L n norm as m and t increase. If the 
boundedness of / is dropped from the assumptions, a linear function would 
be allowed, but assumption (IV) would force it to be decreasing. Since the 
spectrum of —A is positive, (*') is then solved uniquely with u = (which 
is achieved vacuously in the procedure). An example of a function satisfying 
(I)-(IV) is 



/(*) = cor\x) 



whose derivatives are 



f(x) = —±- and f"(x) 2 ''' 



l + x 2 J w (l + x 2 ) 2 ' 

Similarly, if e > 0, A > 0, and h, k e M, then 

Acor 1 ^LzJl^j + k 

represents a family of functions, each of which satify (I)- (IV). A subset of this 
family, given by 

f e (x) = -cot- 1 (-)-l, 
7r e 



is of interest since 



f e (x) -» -H as 
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where H is the Heaviside function and 8 is the Dirac delta function and the 
arrows imply at least pointwise convergence and possibly a more refined limit. 
It is natural to ask whether the Newton-imbedding procedure can be carried out 
in a distributional setting with / = — H and whether f f produces a mcaningul 
approximation to the Heaviside function for small e. More generally if V is the 
class of functions which satisfy (I)-(IV), it is of interest as to which functions 
exist in a suitable closure of V. In this case, 'suitable closure' can be taken to 
mean one whose functions allow for the application of the Newton-imbedding 
procedure in possibly a distributional or more general setting, and produce a 
solution which can be approximated by applying the procedure to a function in 
V. 
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